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Acoustic wave propagation in liquid media containing many parallel air-filled cylinders is con- 

, sidered. A self-consistent method is used to compute rigorously the propagation, incorporating all 

• orders of multiple scattering. It is shown that under proper conditions, multiple scattering leads to 

^\ ' a peculiar phase transition in acoustic propagation. When the phase transition occurs, a collective 

I , behavior of the cylinders appears and the acoustic waves are confined in a region of space in the 

pH ' neighborhood of the transmission source. A novel phase diagram is used to describe such phase 

r-{ I transition. 
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I When propagating through media containing many scatterers, waves will be scattered by each scatterer. The 
^ scattered wave will be again scattered by other scatterers. Such a process will be repeated to estabhsh an infinite 

(~| ■ recursive pattern of multiple scattering, effectively causing the scattering characteristics of the scatterers to change. 
C^i Multiple scattering of waves is responsible for many fascinating phenomena including modulation of ambient 
J-j ' sound at ocean surfaces |^], acoustic scintillation from turbulent flows white paint, random lasers [Q, electrical 
^ I resistivity, and photonic band gaps in periodic structures More interesting, perhaps, under proper conditions 
^ ^ ■ multiple scattering leads to the unusual phenomenon of wave localization, a concept introduced by Anderson Q 
O to explain the conductor-insulator transition induced by disorders in electronic systems. That is, the electrical 
CZJ ■ conductivity can be completely blocked and electrons remain localized in the neighborhood of the initial emission site 
. ^ ' multiple scattering of electronic waves by a sufficient amount of impurities in solids. By analogy, it has been 

c/3 , conjectured that similar localization effect may also exist in the transmission of classical waves in randomly scattering 
media. 

Considerable efforts have been devoted to propagation of classical waves in random media. Localization effects 
have been reported for microwaves in 2D random systems M, for acoustic waves in underwater structures 0], and 
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for light ||9[| . Research also suggests that acoustic localization may be observed in bubbly liquids [g^ 11 1 . Despite the 
CO efforts, however, no deeper insight into localization can be found in the literature, as suggested by Rusek et al. [ p2| . 
^ ^ The general cognition is that enhanced backscattering is a precursor to localization and waves are always localized 
in 2D random systems. Important questions such as how localization occurs and manifests remain unresolved. It is 
also unknown how many types of localization exist. It is believed that localization is a phase transition. However, 



o 
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, how to characterize such a phase transition has not been considered in the literature. A deeper question may concern 
' whether wave localization corresponds to a symmetry breaking and whether the collective behavior often seen in phase 
0^ . transitions such as superconductivity exits. This paper attempts to shed light on these questions. 
0^ . In this paper, we present a rigorous study of acoustic propagation in liquid media containing many air-filled cylinders. 
c/3 ' The approach is based upon the self-consistent theory of multiple scattering [p^ and has been used previously to study 
O , acoustic localization in bubbly liquid media [ pd[ | . Wave propagation is expressed by a set of coupled equations and is 
lAl ' solved rigorously. We report a new phase transition in acoustic localization in 2D random liquids that not only waves 
are confined near the transmitting source but an amazing collective phenomenon emerges. Unlike previous 2D cases, 
waves are not always localized. An essential component of localization in this case arises from the natural resonance 
and the collective behavior of the cylinders. 

We consider a 2D acoustic system in the cylindrical coordinates shown in Fig. 1. Assume that N uniform air- 
cylinders of radius a are placed randomly in parallel in water perpendicular to the x — y plane, forming a random 
cylindrical array. All cylinders are placed within a circle of radius R; the cylinder numerical density is N/{'kB?'). The 
fraction of area occupied by the cylinders per unit area (3 is No?' / B?'] the average distance between nearest neighbours 
is therefore (7r//3)^/^a. No cylinder is located at the origin, where a fine source is placed, and no two cylinders can 
occupy the same spot, i. e. hard sphere approximation. Experimentally, the air-cylinders can be any gas enclosure 
with a thin insignificant elastic layer, like the Albunex bubbles used in echocardiography [|l^. We investigate the 
acoustic propagation in such random liquid media. 
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A unit acoustic line source, transmitting monochromatic pressure waves, is set at the origin, also perpendicular to 
the X — y plane. Due to the large contrast in acoustic impedance between air and water, the air-cylinders are strong 
acoustic scatterers. Unlike the spherical bubble case, the absorption caused by thermal exchange and viscosity is 
unimportant here and can be ignored ||l5[] . 

Multiple scattering in such systems can be computed exactly using the well-known self-consistent method ]l3t and 
the matrix inversion scheme In the scheme, the multiple scattering can be solved in terms of the response function 
of individual scatterers. 

The response function of the single cylinder denoted by 11^ is readily computed by the modal-series expansion in 
the cylindrical coordinates [|l6|,0; for the present case. Hi = 112 = • • • = n^v- It is found that strong response from a 
single cylinder occurs for ka ranging roughly from 0.0005 to 0.5, and a resonant peak is located at ka around 0.005; 
here k is the conventional wavenumber and a is the radius of the cylinders. In this range of ka, to which our attention 
is restricted, the radial pulsation of the cylinders, i. e. the monopole mode, dominates the scattering. In this case, the 
scattering from any single cylinder is isotropic in the x — y plane. The scattering function is plotted against frequency 
in Fig. 2, where both real and imaginary parts of the scattering function are plotted. The scattering cross section 
of a single cylinder reaches a peak. For frequencies above the peak, although it reaches a minimum, the scattering 
function does not vary much. 

Consider that the cylinders are located at [i — 1, 2, .... N) in the cylindrical coordinates. The scattered wave Ps 
from each cylinder, the i-th say, is a linear response to the incident wave composed of the direct incident wave from 
the source and all scattered waves from other cylinders. 

Similar to the 3D case , the scattered wave from the z-th cylinder is therefore 

I 

zTTHl,'\k\f-n\), (1) 

where the second term on the RHS refers to the multiple scattering, and po is the direct incident wave from the source 
and equals the usual 2D Green's function for a unit line source, i. e. iTriJQ^' (fc|r|); here Hq^"^ is the zero-th order 
Hankel function of the first kind. To solve for ps{ri,rj), we set fat one of the scatterers. Then Eq. ([l]) becomes a 
set of closed self-consistent equations and can be solved by matrix inversion. Once ps{ri,rj) is determined, the total 
wave at any space point is given by 

N 

P(r) = Poir) +^Ps{f, n). (2) 

i=l 

First consider wave transmission through the media as a function of frequency. We plot the ratio |p(f)P/|po(^)P in 
the limit |r| — > oo, giving the far field transmitted intensity. A set of numerical computations has been performed for 
various area fractions f3, numbers iV, and dimensionless ka. All parameters are non-dimensionalized. Distances have 
units of the unspecified constant a and functions are dependent only on the dimensionless parameter ka rather than k 
and a independently. Figure 3 presents one of the typical results for the transmitted intensity in a given direction as 
a function of frequency in terms of ka the air-cylinder distribution shown in Fig. 1. It is shown that the transmission 
is significantly reduced from ka = 0.006 to ka — 0.02; little energy is transmitted through. Outside this frequency 
region, however, there is no significant propagation inhibition; transmitted intensity varies around the value one. This 
result holds for any directions of transmission, i.e. it is representative of the frequency dependence of the total power 
escaping to infinity. In line with the 3D situation [0, the regime in which the transmission is blocked indicates 
acoustic localization and is referred to as the localization regime. When N is increased while keeping /3 constant, 
the transmitted intensity for ka between 0.006 and 0.02 will decrease as exp{—{N/ f3y^^ / (La/a)), from which the 
energy localization length La can be estimated in terms of the cylinder radius a. Note from Fig. 3 that at certain 
frequencies, the transmission is enhanced. Such enhancement can be due to strong scattering of individual bubbles, 
multi-refiection at the sample boundary, or the strong multiple scattering among scatterers in certain directions, but 
it is not due to a global collective behavior of the scatterers. We also note that the localization does not occur at the 
resonance peak located at around ka = 0.005. 

To further explore the behavior of acoustic localization in such systems, we investigate the response of each individual 
cylinder to the incident waves. Upon incidence, each cylinder acts effectively as a secondary source. The scattered 
wave can be regarded as the sum of the radiated waves from these secondary sources. From Eq. the contribution 
from a given cylinder may be rewritten as 

p,if,n)^i7rA,H^''\k\f^n\), {i^l,2,...,N), (3) 
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where the complex coefficient Ai denotes the effective strength of the secondary source and is computed incorpo- 
rating all multiple scattering effects. Note that without the scatterers, obviously we will have Ai = IliPo{ri) = 

We express Ai as [^4^1 exp(z6'i): the modulus Ai represent the strengths, whereas di refer to the phase of the effective 
secondary sources. We assign a unit vector, Ui, termed the phase vector hereafter, to each phase 0i, and represent 
these phase vectors by a phase diagram in the x — y plane: The phase vectors are written as Ui = cosOiCx + sin^^ey, 
where and ey are unit vectors in the x and y directions respectively; in the phase diagram, the phase vector Ui is 
located at the cylinder to which the phase 9i is associated. 

Setting the phase of the initial driving source to zero, numerical experiments are carried out to study the behavior 
of the phase vectors and energy spatial distribution as the incidence frequency varies. We observe an amazing pattern 
transition for the phase vectors, with which the wave localization is correlated. The left column of Fig. 3 shows the 
phase diagram for the phase vectors. The right column shows surface plots of the corresponding secondary source 
strength magnitudes |ri||A,;p, giving the energy distribution scaled to eliminate geometrical spreading effects. The 
same cylinder array as in Fig. 1 has been used; three frequencies below, within, and above the localization regime are 
chosen according to the results in Fig. 3. In this particular case, we have set (3 = 10~^ and N = 200. 

We observe that for frequencies below about ka = 0.006 i. e. for frequencies on the left side of the localization 
regime indicated by Fig. 4, there is no ordering in the phase vectors Ui. The phase vectors point to various directions. 
The energy distribution is extended in the x — y plane, and no acoustic localization appears. These are illustrated by 
the case of ka = 0.0044 in Fig. 4. 

As the frequency increases, an ordering in the phase vectors and the energy localization starts to appear. There 
is a transition period in which the phase vectors point to either the positive x-axis or the negative x-axis, as shown 
by the case of ka = 0.008. For frequencies fully within the localization regime indicated by Fig. 3, all phase vectors 
point to the same direction, implying that all cylinders oscillates completely in phase. As indicated by the cases with 
ka = 0.012 and 0.018 in Fig. 4, all phase vectors point to the negative x-direction, exactly out of phase with the 
transmitting source. Such collective behavior allows for efficient cancellation of incoming waves |17| . In this case, 
the wave energy is localized near the transmitting source. It decays about exponentially in all directions, setting the 
localization length to be around 26a, much smaller than the range enclosing the random cylinder array; for /3 = 10"'^ 
and N = 200, the range R is around 450a. The localization behavior is independent of the outer boundary and always 
surfaces for sufffciently large /3 and N. In the field theory, such a collective phenomenon is an indication of a global 
behavior of the system and may imply a symmetry breaking and appearance of a kind of Goldstone boson. 

The non-localization at ka = 0.0044 is not caused by the finite sample size, because no indication of energy 
localization appears as we increase the sample size. If the wave were localized at this frequency, the localization 
length would be shorter than that for ka — 0.012, as the scattering is stronger (Fig. 2). The localization would be 
stronger; this is not observed. 

When the frequency increases further, moving near the upper transition edge, the global in-phase ordering starts 
to disappear. This is illstrated by the case of ka — 0.0216. When moving out of the localization regime, the in-phase 
ordering disappears completely. In the meantime, the wave becomes extended again. This is shown by the example 
ka = 0.04 in Fig. 4. 

In Fig. 4, we also plot the phase portraits for three frequencies {ka — 0.002, 0.0034 and 0.036) at the enhanced 
transmission shown by Fig. 3. It is shown that there is no apparent phase ordering in these cases. The figure also 
show that the enhancement in one direction does not necessarily imply enhancement in all directions. 

To quantify the phase transitions shown in Fig. 4, we may define two tentative order parameters. For the in-phase 
behavior, the magnitude of the summation of all phase vectors normalized by N may be regarded as a phase order 
parameter, called order parameter 1. For the mixed ordering shown in the transition regime, we may define the order 
parameter as {\ux{i)\ — \uy{i)\), termed order parameter 2. These order parameters are plotted as a function of 

frequency in Fig. 5. When all cylinders oscillate in phase, the order parameter 1 is close to one, while in the transition 
regime, the order parameter 2 may be more proper in describing the phase transition. 

The behavior of the system is a result of the interplay between the single cylinder scatter function and the many 
body interactions. Without multiple scattering, the localization and coherence phenomenon does not occur. Both 
the phase and strength of the scattering function determine how a multiply scattered wave propagates through the 
medium. Exactly how the interplay between the scattering function and the multiple scattering gives rise to the 
localization and the global coherence behavior remains an open question and may only be understood by an analytic 
analysis which is under our consideration. 

The localization or the transmission drop starts at a frequency slightly above the resonant peak, about the minimum 
of the real scattering function. Throughout the localization regime, the scattering function varies mildly. The 
localization or the transmission drop disappears at some frequency which does not appear to be dependent on the 
peak or the dip of the scattering cross section. The phase transition at the lower frequency edge seems to be related 



3 



with the small dip in the scattering cross section or in the imaginary part of scattering function, while the transition 
at the higher end to the decreasing wavelength. Once the cylinders become too many wavelengths apart, despite 
a slight increase in the strength of the scattering cross section, phase coherence can no longer be maintained and 
localization is destroyed. 

We have modified the scattering function manually, and observed that transmission is sensitive to it. If, however, 
the scattering function is held constant at its value for a frequency within the localized range, and the transmission 
calculated versus frequency using this constant value, the phase transitions and a similar transmission behavior are 
still observed. 

The following considerations support and elucidate the localization behaviour and phase transition observed. (1) 
The localization is not due to dissipation. Acoustic absorption is negligible in the cases considered, and no mechanism 
for absorption has been included in our model. (2) By varying the cylinder number N while keeping the area fraction /3 
constant, it can be shown that the localization or non-localization behavior is qualitatively unchanged, thus not caused 
by the boundary of the cylinder arrays. (3) The above localization phenomena are caused by multiple scattering. 
When we manually turn off the multiple scattering from Eq. (|l|), the localization disappears. (4) The localization 
behavior also disappears when the cylinder concentration reaches a sufficiently low level. The localization range in 
Fig. 3 is narrowed, tending to vanish. (5) When the wave is not localized, the energy distribution varies as the 
random placement of the cylinders changes at a given concentration (3. Once localized, the localization behavior will 
not change as the cylinder placement varies. 

In summary, we have demonstrated a new phase transition in acoustic propagation in 2D random arrays of air- 
cylinders in water. When the concentration of the air-cylinders exceeds a certain value, acoustic waves become 
localized near the transmitting source within a range of frequencies. The results indicate that the wave localization is 
related to the collective behaviour of the air-cylinders in the presence of multiple scattering. We take the view that 
the coherent behavior allows for efficient cancellation of the source wave, giving rise to localization, but localization 
and coherence occur hand-in-hand and the statement that the localization is due to the coherence or visa-versa may 
be just a matter of perspective. Although these properties may not hold in general, the fact that they do for resonant 
air-cylinders makes these scatterers ideal for theoretical and experimental localization studies. 

The work received support from National Science Council of ROC and from National Central University in the 
form of a special scholarship to EH. 
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FIGURE CAPTIONS 



Fig. 1 Top view of a random distribution of air- cylinders. The circles refer to the air-cylinders (not to scale). A 
straight line source is placed at the origin. 

Fig. 2 The scattering function versus frequency. 

Fig. 3 Transmission versus frequency in terms of ka. 

Fig. 4 Left column: Phase diagram for the two-dimensional phase vector defined in the text. Right column: Spatial 
distribution of acoustic energy (arbitrary scale). 

Fig. 5 'Order parameters' versus frequency. 
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